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ABSTRACT 

We have used multi-frequency follow-up observations of a sample of extragalactic sources 
from the 9C survey at 15 GHz to make deductions about the expected source population at 
higher radio frequencies, such as those in the lower frequency bands of the Planck Surveyor 
satellite. In particular, we have made empirical estimates of the source counts at 22 GHz, 
30 GHz, 43 GHz and 70 GHz and compared these with both known data and current theo- 
retical predictions. We have also made an estimate of the count at the ALMA frequency of 
90 GHz, with a view to assessing the possible population of point sources available for the 
phase calibration of that instrument. 
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1 INTRODUCTION 

A major problem for many cosmic microwave background (CMB) 
experiments, such as the Planck Surveyor mission, is the estima- 
' tion of the contaminating effect of extragalactic radio sources on 
the measurement of the CMB anisotropies. In order to model the 
expected response from a CMB instrument, with a view to develop- 
ing appropriate analysis techniques, it is essential to have a realistic 
assessment of the source count at the relevant frequency. There has 
1 been important work on modelling the source counts at high fre- 
■ quencies, such as that of Toffolatti et al. (1998) and De Zotti et al. 
(2005), but few direct observational measurements are available. 
In particular, in the bands of the Planck Low Frequency Instru- 
ment (LFI) - i.e. 30 GHz, 44 GHz, 70 GHz - there is, at present, a 
lack of data from appropriate wide-field surveys, although several 
such surveys are under way: AT20G (Ricci et al. 2004, Sadler et al. 
2006), OCRA (Browne et al. 2000), RATAN-600 (Parijskij 2004). 
The 9C survey (Waldram et al. 2003) at the comparatively close 
frequency of 15 GHz is therefore very important for making some 
predictions. 

In this paper we have made use of data from the 'simul- 
taneous' multi-frequency 9C follow-up observations described in 
Bolton et al. 2004 (paper 1). There have already been two other 
papers based on these data: Bolton et al. 2006a (paper 2), which 
presents 5-GHz MERLIN and VLBA observations of a set of 36 
compact sources, and Bolton et al. 2006b (paper 3), which reports 
results from a study of the 15-GHz variability of 51 sources. Here 
we have taken a complete sample of 1 10 sources above 25 miy at 
15 GHz and used the spectra over the range 1.4 GHz to 43 GHz to 
make some deductions about the expected source population in the 
Planck LFI bands. Our approach is entirely empirical and there is 
no attempt to model the sources. 

We have also used our data to make a prediction of the count at 



90 GHz in order to estimate the number of point sources available 
as possible phase calibrators for the ALMA telescope. 

In section 2 we describe the sample and include some discus- 
sion of variability. Section 3 explains the principles of our method 
and the assumptions made. In section 4 we examine the spectral 
index distributions over the range 1.4 to 43 GHz and in section 5 
we discuss the interpolation and extrapolation of the source spectra 
to 30, 70 and 90 GHz. Section 6 presents the results of our cal- 
culations and section 7 attempts to estimate the errors involved. 
In section 8 we test our method by comparing our predictions for 
1.4 GHz and 4.8 GHz with the known counts at these frequen- 
cies. Section 9 presents estimates of the source counts at 22 GHz, 
30 GHz, 43 GHz, 70 GHz and 90 GHz and compares them with 
those from the available high frequency data and with the theoreti- 
cal models of De Zotti et al. (2005). Section 10 refers to the 90 GHz 
result and its implication for ALMA phase calibration. Finally, sec- 
tion 11 is a general discussion of the reliability of our method and 
section 12 is a summary of our conclusions. 

We note here that our definition of spectral index, a for flux 
density, 5, and frequency, v, is 5 = v^". 



2 THE SAMPLE 

We assembled a flux-limited sample of 121 sources, complete to 
25 mjy at 15 GHz, from three areas of the 9C survey. All the 
sources had a complete set of simultaneously measured flux den- 
sities at frequencies of 1.4 GHz, 4.8 GHz, 15.2 GHz, 22 GHz and 
43 GHz apart from 11 of them, for which some observations were 
missing. Since the reasons for the omissions were totally unre- 
lated to the source characteristics - i.e. they were due to weather 
or scheduling problems - we have simply omitted these sources 
from the sample and calculated the final effective area from a fit of 
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Table 1. The 1 10 sources in our sample with their simultaneous flux densi- 
ties at 1.4, 4.8, 15.2, 22 and 43 GHz: S1.4, S4.8, 5i5.2, S22, S43 in mjy 
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Table 1 - continued 
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the known 15 GHz source count (see below). We thus have a 'snap- 
shot' of sets of instantaneous spectra for a complete sample of 1 10 
sources, in an area of approximately 130 deg^, as hsted in Table 1. 
(For more details of individual sources see papers 1 and 2.) 

The extent of variability at 15 GHz, as illustrated by our sam- 
ple, can be seen in Figure 1 which shows how the sources have 
changed in flux density beween the time of the original observa- 
tions and the time of the follow-up observations. The sources in 
the original selection had observation dates ranging from Novem- 
ber 1999 to September 2001 and the follow-up observations had 
dates in January, November, December 2001 and January, May 
2002. The maximum time interval between the original and follow- 
up observations of any source was 30 months and the minimum 
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above 100 mjy. The implications of this are discussed in later sec- 
tions of the paper. 
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Figure 1. Plot of 15 GHz follow-up flux densities versus original flux den- 
sities, showing the scatter about the line of unit slope. The dotted lines cor- 
respond to a flux density of 25 mJy on each axis. 



was 2 months. The source showing maximum variability was 
J1514-I-3650 which had risen from 40 mJy to 95 mJy in a period 
of 19 months. Since the time intervals between our original and 
follow-up observations differ widely, Figure 1 can give only a gen- 
eral indication of the extent of variability at 15 GHz. (For more 
detailed work see paper 3.) 

Variability presents a problem for the selection of a flux lim- 
ited sample, since close to the lower limit there will be preferential 
selection of those sources which were above their mean flux den- 
sity values at the time of observation rather than below. Hence, for 
our sample. Figure 1 shows an excess of sources with flux densi- 
ties below 25 mJy in the follow-up observations, and this is, in fact, 
unavoidable. 

Another feature of our sample is that it contains significantly 
fewer sources above 100 mJy than are expected from the known 
15 GHz source count (Waldram et al. 2003). This became apparent 
when we were calculating the effective sample area from a fit of the 
known count to the flux density distribution in the original sample. 
We found that in the range from 100 mJy to the maximum value 
of 665 mJy there were only 5 sources and we therefore used only 
those below 100 mJy in our fit. (In fact, a value for this area is 
not required for our predictions and is used only in Section 7 and 
Figure 7 for calculating the incomplete counts derived directly from 
the sample, by way of comparison.) 

Similarly, in the follow-up sample there are only 6 sources 
above 100 mJy, the highest value being 540 mJy. Taking account of 
the uncertainty in the area and the Poisson errors, we estimate that 
we should expect at least twice this number. The short-fall may be 
be due to the fact that, in the original selection of the 9C fields 
for the purpose of CMB observations, there was some bias against 
regions predicted to contain very bright sources (see Waldram et al. 
2003). However, our area here is quite small and, ideally, a much 
larger area is required to sample adequately the population above 
100 mJy. 

Thus, although for the purpose of this work we can assume 
that our original sample represents the typical characteristics of a 
complete sample of sources in the flux density range 25 mJy to 
100 mJy, there is less certainty in the under-sampled upper range. 



3 EMPIRICAL ESTIMATION OF THE SOURCE 
COUNTS 

We now consider whether, knowing the source count at 15 GHz, 
it is possible to estimate the source counts at higher frequencies 
from these data. It is not possible to do so directly because a com- 
plete sample would be far too small. For example, if we assume an 
extreme rising spectral index between 15 and 43 GHz, Ujj, of —1, 
our sample at 43 GHz is complete to only ~ 70 mJy, providing only 
9 sources. It is however possible to use the 15 GHz count and our 
distribution to estimate the count at 43 GHz, if we make certain 
assumptions about the ajj distribution in our flux-density range. 

Consider first a source population such that each source has 
the same spectral index a between the two frequencies Vi and V2, 
or Svi = rSv2 where r = (V1/V2) and let us assume that the 
differential count at Vi has the form 

nvi (5) = diVvi /dS = Ay^ S'^ 

where Av, and b are constants. 

In order to find the differential count at V2 we consider the cor- 
responding integrated counts, A'vi(> S) and iVv2(> '^)- fr^" 
quency Vi (for b 7^ 1) 

and so at frequency V2 

Wv, (> S) = Nv, (> rS) = Av, ( 1 - b)-l (r5) 

This means that the differential count at V2 becomes 

(S) = dNvJdS = S'^ . 

In practice, however, we know that the source population 
spans a range of spectral indices and we now make the following 
assumptions: first, that our sample provides a tjfpical distribution of 
spectral indices and secondly, that this distribution is independent 
of flux density. In our sample of m sources we know the spectral 
index a, and the corresponding value of r, r,-, for each source, and 
so can calculate ki = r?"**, and hence 

where 

1 m 1 m 

i—i i—i 

(See Condon 1984 and Kellermann 1964 for similar analyses.) 

We see that these assumptions lead to a form of the count at 
V2 with the same exponent as the count at Vi but with a different 
pref actor. We can define an effective value for r, r^, such that 

where re would be the value of r for all sources, if they all had the 
same spectral index. Using we can then estimate the flux density 
range over which the count derived for frequency V2, from a knowl- 
edge of the covmt at frequency Vi (i.e. at 15 GHz), may be assumed 
to be valid. For this purpose wc define values at V2: S^jn = '^min/'^e, 
Smax =S^/re md Sc = S^^ / r^,, where and^^^j^ are the min- 
imum and maximum values in the original 15 GHz sample, 25mJy 
and 665 mJy respectively, and S^^ is the upper 'completeness' value 
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Figure 4. Plot of versus afj, 
0(^2 = «i5- Sources with 15 GHz (follow-up) flux densities above 100 mjy 
are ringed. 



Table 2. Table showing the distribution of numbers of sources with respect 
to the median spectral index (ttmed) and the median 15 GHz follow-up flux 
density (Smed = 43.0mjy). See Figure 2. 
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of 100 mJy above which we know the data are sparse. We assume 
that, although our predictions may be reliable in the range Smm to 
5c, they will be less so in the range Sc to 5niax- 



4 SPECTRAL INDICES IN RANGE 1.4 TO 43 GHZ 

Figure 2 shows histograms of the spectral index distributions in the 
range 1 .4 to 43 GHz. In the case of the 43 GHz observations there 
were 7 non-detections of sources. For these we have set the flux 
density value equal to the noise level and the corresponding spec- 
tral indices are shown shaded in the figure. It can be seen that the 
median spectral index increases with increasing frequency, being 
0.37 for a}^^, 0.47 for al\, 0.54 for a^^ and 0.89 for af^. In Fig- 
ure 3, in order to investigate the dependence of spectral index on 
flux density, we have drawn scatter plots of the 15 GHz follow- 
up flux density versus spectral index for the four frequencies and 
added lines showing the median values of the two quantities. There 
is no apparent correlation visible and this is confirmed by Table 2 
which gives the number of sources in each of the four quadrants of 
each plot. We conclude that, considering the Poisson errors in these 
small numbers, the differences are not significant, although, as we 
have seen in Section 2, the distributions are undersampled above 
100 mJy. 

In Figure 4 we have drawn a scatter plot of versus a^y 
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Figure 5. Spectrum of J1530+3758, showing the known data points and the 
points calculated from a quadratic fit of logS with log v. 



Sources with 15 GHz (follow-up) flux densities above 100 mJy 
are ringed. It can be seen that only three sources have rising spec- 
tra from 22 to 43 GHz: these are JOOlO-l-2854, J0019-I-2817 and 
J1520-I-4211 and have 15 GHz (follow-up) flux densities of 69.2, 
17.4 and 53.8 mJy respectively, all, in fact, less than 100 mJy. 

For further discussion of the reliabiUty of the spectral index 
measurements see section 11. 



5 INTERPOLATION AND EXTRAPOLATION OF THE 
SOURCE SPECTRA 

In order to estimate the value of K for the frequencies 30, 70 and 
90 GHz, we have needed to interpolate and extrapolate the spectra 
of the individual sources by an appropriate fit to the existing data. 
From the 'colour-colour' plot in Figure 4 we see that, out of a total 
of 1 10 sources, 91 have both ajg and > and we assume that 
these are steep spectrum sources whose spectra have already turned 
over at frequencies below 43 GHz. Of these, 75 have ttj^ > 01^5 and 
so any fit should attempt to take account of this steepening of their 
spectra (see, for example, the plot for J1530-I-3758 in figure 5). 

We have tried three types of fit of log 5 with logv : (a) a 
quadratic fit for all sources using the values at 15.2, 22 and 43 GHz 
, (b) a similar quadratic fit for the 75 sources mentioned above but 
a Unear fit for the remainder, using the 22 and 43 GHz values and 
(c) a linear fit for all sources, using the 22 and 43 GHz values. The 
linear fit is equivalent to the assimiption that a™ and are both 
equal to a^2- 

In each case, our procedure was to apply the appropriate fit 
to the known spectral points of each source and calculate the flux 
densities at 30 GHz, 70 GHz and 90 GHz by interpolation or ex- 
trapolation. We could then find the values of r, at each frequency 
and hence the corresponding values of K. 

We examined the results of the three types of fit for a selection 
of sources from the different areas in the 'colour-colour' plot. Fit 
(a) catered well for the 75 sources with steepening falling spectra 
but produced some obvious anomalies for a number of the flatter 
spectra. Fit (b) appeared to produce much more appropriate results 
for all types of spectra, while fit (c) did not take sufficient accoimt 
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spectral index {1 5 to 22 GHz) 



spectral index {1 5 to 43 GHz) 



Figure 2. Distributions of spectral index: (x{ 4, a^, afj, a\l. The shaded area indicates the 7 non-detections at 43 GHz where the flux densities have been set 
equal to the noise level. 



Table 3. Table showing the values of K and their uncertainties for the four 
frequencies 1.4, 4.8, 22 and 43 GHz. 



Frequency 


K 


K 


K 


GHz 


all sources 


Si5 < 43 mjy 


Si5 > 43 mJy 


1.4 


5.38 ±0.53 


5.51 ±0.77 


5.25 ±0.73 


4.8 


2.15±0.13 


2.23 ±0.18 


2.08 ±0.18 


22 


0.80±0.02 


0.78 ±0.03 


0.82 ±0.03 


43 


0.44 ±0.04 


0.42 ±0.05 


0.45 ±0.05 



Table 4. Table showing the values of K and their uncertainties for the three 
frequencies 30, 70 and 90 GHz. See section 5. 



of the steepening falling spectra. It was therefore decided to use the 
values of K from fit (b), but the results from all three types of fit are 
included in Table 5. 

We should emphasise that this method of interpolation and 
extrapolation is purely empirical and no attempt has been made to 
model the sources or to examine every source individually. 







K 

all soitfces 


K 

Sis < 43 mJy 


Sis > 43 mJy 


30 GHz 












(a) 


0.62 ±0.03 


0.60 ±0.04 


0.65 ±0.05 




(b) 


0.61 ±0.03 


0.59 ±0.04 


0.63 ±0.04 




(c) 


0.59 ±0.03 


0.57 ±0.04 


0.61 ±0.04 


70 GHz 












(a) 


0.25 ±0.03 


0.25 ± 0.04 


0.25 ± 0.05 




(b) 


0.29 ±0.04 


0.29 ± 0.06 


0.30 ±0.07 




(c) 


0.33 ±0.04 


0.31 ±0.05 


0.34 ±0.07 


90 GHz 












(a) 


0.19 ±0.03 


0.21 ±0.04 


0.18 ±0.04 




(b) 


0.25 ±0.05 


0.25 ±0.06 


0.26 ±0.08 




(c) 


0.29 ±0.05 


0.27 ±0.06 


0.31 ±0.08 
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Figure 3. Scatter plots of the 15 GHz foUow-up flux density versus spectral index for: a[ 4, a4^8, afj, a^j. In each case the dashed lines show the median flux 
density and median spectral index. See Table 2. 




Table 5. Table showing the results of our calculations for the frequencies 
from 1.4 GHz to 90 GHz 



Frequency 


K 


A 


re 


5min 




5c 


GHz 




Jy-'sr-i 




mjy 


mJy 


mJy 


1.4 


5.38 


274 ±30 


0.23 


110 


2900 


435 


4.8 


2.15 


110±9 


0.51 


49 


1300 


195 


15.2 


1.00 


51±3 


1.00 


25 


665 


100 


22.0 


0.80 


41±2 


1.21 


21 


550 


85 


30.0 


0.61 


31±2 


1.54 


16 


430 


65 


43.0 


0.44 


22±2 


2.04 


12 


330 


50 


70.0 


0.29 


15±2 


2.93 


9 


230 


35 


90.0 


0.25 


13±3 


3.34 


7 


200 


30 



Figure 6. Distribution of ki values for 43 GHz. 



6 THE CALCULATIONS 

For each frequency we have calculated the value of ki for each 
source, where ki = r\^^, and then taken the unweighted mean of 
the distribution, K = A:, . This was repeated with the sam- 

ple divided into two groups, one with 5] 5 < 43 mJy and the other 
with 5i5 > 43 mJy, 43 mJy being the median flux density of the 
15 GHz follow-up observations. The results for the frequencies 1.4, 
4.8, 22 and 43 GHz are shown in Table 3. In Table 4, for 30, 70 and 



Some estimates of the source counts at Planck Surveyor frequencies, using the 9C survey at 15 GHz 1 



90 GHz, we have also included the results from the three meth- 
ods of spectral interpolation/extrapolation, (a), (b), (c), described 
above. In each case the uncertainty quoted is the error in the mean, 
or o/x/w , where ct is the standard deviation of the fc, distribution 
and m is the total number of sources. The distributions are neces- 
sarily skewed; an example is shown in Figure 6, for 43 GHz. 

Table 5 shows the results of our calculations. We take the 
coimt at 15 GHz to be 

(Waldram et al. 2003). At another frequency v the exponent 
b ( = 2.15) remains the same but the prefactor A becomes A" x 5 1 . 
For each frequency we have calculated r^, and also the values of 



7 ERROR ESTIMATES 

There are a number of factors which contribute to the imcertainty 
in our values in K and A. As well as the error in the mean of the dis- 
tribution ki = r]^^ , (o/i/m), there is also the error in the original 
15 GHz source count: i.e. for b = 2.15 ± 0.06 the error in the pref- 
actor (A at 15 GHz) at the centre of the data is ~ 5%. Combining 
these gives the uncertainties in A shown in Table 5. 

These are the minimum estimates for the uncertainties. The er- 
rors hardest to quantify in our source count predictions lie in the 
assumption that the spectral index distribution is independent of 
flux density. We have seen from Table 2 and Figure 3 that we can- 
not detect a dependence within our current sample. Similarly, in 
Tables 3 and 4, the differences in the K values for sources with 
5i5 < 43 mjy and those with S\i > 43 mjy are not significant 
given the intrinsic errors. However, with only 6 sources above 100 
mJy, our distributions of ki may be skewed towards values of K 
more appropriate to the lower rather than the higher flux densi- 
ties. Also, it is possible for sources outside our selected flux den- 
sity range at 15 GHz, with a different spectral index distribution, to 
contribute to the predicted counts at another frequency, even within 
our estimated range of Si-^m to 5max for that frequency (see section 
11 for further discussion). 

We need to investigate how far the assumption of the indepen- 
dence of spectral index and flux density is a useful approximation 
for our present work. One way of testing our procedure has been 
to apply it to our data for 1.4 and 4.8 GHz, since we already have 
measurements of the source counts at these frequencies from earUer 
surveys (see section 8). 

We can also gain some insight into the reliability of our 
method by comparing our predicted counts for 22 and 43 GHz with 
the incomplete counts derived directly from the sample, as in Figure 
7. Here we show the two types of error on each data point: one is a 
systematic error due to the uncertainty in deriving the sample area 
from a fit to the known 15 GHz count (see section 2) and the other 
is the usual random Poisson error. We see that the predicted counts 
are close to the direct counts at the higher flux densities whereas 
at the lower flux densities there is a marked fall-off in the direct 
counts. This is as expected, since the original sample was complete 
to only 25 mJy. We have calculated the values of S above which 
we might expect the direct counts to be complete, assuming an ex- 
treme rising spectral index, af^ or a^j, of —1. These are 36 mJy 
for 22 GHz and 71 mJy for 43 GHz, and we see they are consistent 
with these plots. 
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Figure 8. Prediction for 1 .4 GHz count, where the two dashed lines show 
the uncertainty in A from Table 2. The vertical dotted lines indicate Smm 
andSmax 
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Figure 9. Prediction for 4.8 GHz count, where the two dashed lines show 
the uncertainty in A from Table 2. The vertical dotted lines indicate Smin 
andSmax 

8 LOWER FREQUENCY COUNTS 

Figure 8 shows a plot comparing our predicted count at 1.4 GHz 
with a count from the NVSS survey (Condon et al. 1998). We have 
used 1.4 GHz data from the areas of our original 9C survey fields, 
as in Waldram et al. 2003. It can be seen that our prediction lies 
significantly below the data points in the range of flux density over 
which we might expect it to be valid. We estimate that an approxi- 
mate fit to the data, keeping the count exponent at —2.15, is given 
by a value of A of ~ 480, as compared with the predicted value of 
274 ± 30. This is not surprising because we already know that the 
assumption of the independence of spectral index and flux density 
does not hold over wide ranges of frequency and flux density, as 
can be seen, for example, in the 1.4 to 15 GHz spectral index dis- 
tributions in Waldram et al. 2003 and Waldram & Pooley 2004. In 
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Figure 7. The predicted counts for 22 and 43 GHz compared with the incomplete counts derived directly from the sample. The two dashed lines correspond 
to the uncertainty in A from Table 2. The pairs of data points indicate the uncertainty in the area (111 to 145 deg^) and the error bars are the Poisson errors. 
For clarity we have artificially separated the points within the pairs so that they are slightly above and below the corresponding value of S. The vertical dotted 
lines show the values of S above wiiich the direct counts are expected to be complete. 



the latter paper we show that the percentage of inverted spectrum 
sources, with a\^^ < 0, increases with increasing flux density: we 
find that for three samples - i.e. 5 to 25 mjy, 25 to 100 mJy, 100 mjy 
and above - the percentages are 10, 20 and 33 respectively. In our 
current sample the percentage is also 20. This would suggest that 
below the completeness limit of our sample there is a source popu- 
lation with a higher proportion of steep spectrum sources and these 
are consequently contributing an extra component to the count at 
1.4 GHz. 

We have also compared our 4.8 GHz prediction with the count 
from the Green Bank survey, using the data in Gregory et al. 1996, 
as illustrated in Figure 9. This corresponds to a frequency ratio of 
only ^ 3, rather than ~ 10, and we can see that there is agreement 
within the errors over a flux density range of approximately 0.175 
to 1.25 Jy. At lower flux densities our predicted count is too low 
and at higher flux densities too high, which is consistent with the 
trend in spectral index found for 1.4 GHz. 



9 HIGHER FREQUENCY COUNTS 

At frequencies above 15 GHz the only data available with which to 
compare the predictions are those at 31 GHz from DASI (Kovac et 
al. 2002) and CBI (Mason et al. 2003) and at 33 GHz from the VS A 
(Cleary et al. 2005). Other source count data from ATCA (Ricci et 
al. 2004) or WMAP (Bennett et al, 2003), for example, lie outside 
the relevant flux density range. Figure 10 shows the comparison of 
our predicted count with the measured counts from DASI and CBI 
and with data from the VSA. 

From DASI, over the range 0. 1 to 10 Jy, we have taken 
dAf/d53i = (32±7)(53i/Jy)-2-'5±°-2°Jy-isr-i 
and from CBI, over the range 0.005 to 0.05 Jy, 
iV(>53i) = (2.8±0.7)(53i/10mJy)-i0deg-2 
or 

dAf/dS3i =(92± 23)(53i/Jy)-2'»Jy-isr-l . 

The VSA count, over the range 0.02 to 0.1 14 Jy, is fitted by 

dN/dS33 = (21±4.5)(533/Jy)-^-^4jy-^sr-i 
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Figure 10. Predicted 30 GHz count with VSA data and the measured CBI 
and DASI counts. The two dashed lines show the uncertainty in A from 
Table 2. The vertical dotted lines indicate Sniin and 5max. 



We see from Figure 10 that for 5 > 0. 1 Jy our prediction is consis- 
tent with the DASI count and for S < 0. 1 Jy it is consistent with the 
VSA data, but it lies somewhat below the CBI count in the range 

S < 0.05 Jy 

Although no further direct data are available, we can compare 
our measured 9C count and our empirically predicted higher fre- 
quency counts with those from the models of de Zotti et al. (2005); 
these plots are shown in Figure 1 1 . Here the de Zotti models repre- 
sent the sum of the contributions from the three main extragalactic 
source populations: FSRQs (flat spectrum radio quasars), BL Lacs 
(BL Lacertae type objects) and steep-spectrum radio sources. The 
contributions from other types of extragalactic source population 
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Figure 11. Comparisons of our predictions, and also the known 9C 15 GHz count, with the de Zotti models. The two dashed lines show the uncertainty in 
A from Table 2. For 70 and 90 GHz these correspond to extrapolation of type (b), but the error limits for extrapolation of types (a) and (c) are also included 
(shown dotted). The vertical dotted lines indicate Smin and Smax and the vertical arrow marks the value of Sc. 



10 



are assumed to be negligible over the relevant ranges in flux den- 
sity. We see that there is good agreement between our counts and 
the de Zotti models over the appropriate Smin to 5max range for the 
frequencies 15, 20 and 30 GHz, but for 40, 70 and 90 GHz, al- 
though there is good agreement over the lower part of the range, 
below our upper 'completeness' value ^c, there is an increasing di- 
vergence at the higher flux densities, the models lying significantly 
above our counts. 

In Figure 12 we have repeated the comparison of our 43 GHz 
prediction with the de Zotti model, but here we have included the 
separate model components from FSRQs, BL Lacs and steep spec- 
trum sources. Taking appropriate tabulated values from the de Zotti 
counts, we have made the following calculations. At a flux den- 
sity of 316 mjy, a value close to the Smux of 330 mjy, our count 
is only 65 ± 6 % of the de Zotti 40 GHz total count. At this point 
the de Zotti model is dominated by the contribution from FSRQs, 
which amounts to 83 % of the total. It is apparent that our empiri- 
cal approach is predicting significantly fewer flat spectrum sources 
than the de Zotti model in this higher flux density range. 

However, as we have seen, our 15 GHz sample contains few 
sources above 100 mJy and so our predicted count at another fre- 
quency is less reliable in the range 5c to Smux- It is possible that 
with increasing flux density there is a significant shift in the spec- 
tral index distribution towards flatter spectra, even if such a trend is 
not detectable in our data. 



10 PHASE CALIBRATORS FOR ALMA AT 90 GHZ 

One strategy proposed for the ALMA phase calibration is to use 

point sources with flux densities above 20 mJy at 90 GHz and to ex- 
trapolate the phase solutions up to the appropriate target frequency 
(see, for example, Holdaway & Owen, 2005). It is therefore impor- 
tant to investigate whether there will be a sufficient density of such 
sources available at this frequency. Holdaway & Owen have devel- 
oped a simple parametrized model of the source population, and, 
using observed 8.4 GHz and 90 GHz fluxes, estimated the source 
counts as a function of frequency. They have estimated a count 
of about 1800 point sources per steradian brighter than 20 mJy at 
90 GHz. 

As far as our own 90 GHz prediction is concerned, our esti- 
mate of the number of sources in the range 20 to 200 mJy (where 
200 mJy is the value of 5max) is 940 ± 220 sr" ' . However, not aU of 
these sources wiU be of a sufficiently smaU angular size for use as 
calibrators. At 15 GHz we should expect about half the sources to 
be less than 0.1 arcsec in angular diameter (see paper 2), though at 
90 GHz it is likely to be a somewhat higher fraction. Our estimate 
is thus equivalent to approximately one such source in every 6 or 
7 square degrees in this flux density range. 

If this is correct (though it may well be an underestimate), it 
suggests that suitable phase calibrator sources for ALMA may be 
as much as twice as far from the target source as has been assvmied 
hitherto, meaning that either longer slews or longer on-source inte- 
grations will be required to achieve good phase solutions. However, 
given the very high sensitivity of ALMA, this is not anticipated to 
affect significantly the observing efficiency. 



11 DISCUSSION 

As we have seen, our empirical method of predicting the source 
counts depends on two main assimiptions: first, that our sample 



provides a typical distribution of spectral indices and secondly, that 
this distribution is independent of flux density. 

We consider first the reliability of our spectral index measure- 
ments. Since sources at these radio frequencies can be extremely 
variable (see paper 3) it was essential to make our follow-up ob- 
servations simultaneously and the fact that we were able to do so 
is an important element in our method. However, it has meant that, 
at any one time, the VLA measurements were made with the same 
configuration of the telescope for all frequencies, leading to a wide 
variation in the size of the synthesized beam across the frequency 
range (see paper 1). Although we used integrated flux densities at 
all frequencies, it is possible that some flux of the more extended 
sources has been 'resolved out' at the higher frequencies, resulting 
in a bias towards steeper spectra. 

To investigate the effect this might have on our predictions, we 
divided the sample into two groups, those sources which appeared 
point-like at all frequencies (73) and those which were resolved at 
one or more frequencies (37). For each of the frequencies 22 and 
43 GHz, we further divided the groups into 'FS' or flat spectrum 
sources (ajj < 0.5) and 'SS' or steep spectrum sources (a^j > 0.5). 
Our results are shown in Table 6. We see that at both frequencies the 
mean value ki for the resolved SS sources is lower than the value 
for the point SS sources and that at 43 GHz the difference appears 
to be marginally significant. Of course, we cannot tell whether this 
means that the spectra of the resolved SS sources are genuinely 
steeper, which is quite possible, or whether we are actually miss- 
ing flux from over-resolution. We can, however, assign the mean ki 
value of the point SS sources to each of the resolved SS sources, 
and calculate a test value of K for the whole sample for compari- 
son with the K value used in our predictions. We find (see Table 6) 
that these do not differ significantly, meaning that, even if we were 
losing flux of some of the extended sources, the error is likely to be 
small. 

This conclusion is corroborated by the fact that the predicted 
counts at 30 GHz, calculated by interpolating the source spectra 
between 22 and 43 GHz, are in good agreement with experiment. 

We have seen from section 3 that our second main assumption, 
that the spectral index distribution is independent of flux density, 
leads to a form of the predicted count with the same exponent b as 
at 15 GHz but a different pre-factor A. Any variation in the spectral 
index distribution with flux density can be envisaged as resulting 
in a dependence of A on flux density, equivalent to curvature in the 
logarithmic count. 

We have already emphasized earlier in the paper that our data 
are sparse at the higher flux densities. Thus, although we have 
shown that, within our sample, spectral index and flux density ap- 
pear to be independent, we cannot assume that this continues to 
hold above a 15 GHz flux density of 100 mJy. We have, though, 
been able to assemble from a wider area a small un-biassed sample 
of 16 sources (including 5 from the original sample) with 15 GHz 
flux densities in the range 102 to 784 mJy and simultaneous follow- 
up observations at 15, 22 and 43 GHz. We find that they all have 
spectral indices a'^ in the range 0.1 to 1.1, indicating a possible 
flattening of the spectra compared with the values shown in Fig- 
ure 4. A calculation of the value of A at 43 GHz from these 16 
sources gives A = 27 ± 2, as compared with 22 ± 2 in Table 5, 
corresponding to a count at 316 mJy of 79 ± 7 % of the de Zotti 
40 GHz total count, as compared with the 65 ± 6 % quoted in sec- 
tion 9. There is thus some evidence for an increase in the value of 
A and consequently for closer agreement with the de Zotti predic- 
tion; however, given the error estimates, it is scarcely conclusive. 
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Figure 12. Comparison of our 43 GHz prediction with the de Zotti model, showing the separate model components from FSRQs, BL Lacs and steep spectrum 
sources. The two dashed lines show the uncertainty in A from Table 2. The vertical dotted lines indicate Smin and Smax and the vertical arrow marks the value 
of Sc. 



Table 6. Table showing the results of our tests for the effect of possible bias in the spectral index distributions. 'FS' denotes flat spectrum sources (a^j < 0.5) 
and 'SS' denotes steep spectrum sources (Hjj > 0.5). 



Frequency 




point sources 




resolved sources 


K 


K 


GHz 


no.ofFS 


no. of SS mean fcj for SS 


no.ofFS 


no. of SS mean kj for SS 


test value 


value used 


22.0 


40 


33 0.66 ±0.02 


10 


27 0.62 ±0.03 


0.81 ±0.02 


0.80 ±0.02 


43.0 


27 


46 0.30±0.02 


4 


33 0.23 ±0.02 


0.46 ±0.03 


0.44 ±0.04 



It is clear that more data at these higher flux densities would be 
required to detect any significant change in the slope of the count. 



12 CONCLUSIONS 

We have shown that it is possible to use our multi-frequency follow- 
up observations of a sample of sources from the 9C survey at 
15 GHz to make some empirical estimates of the source counts 
at higher radio frequencies. These predictions are important, in 
spite of the necessary limitations of our method, since at present 
there are few direct observations at these frequencies. Our data, al- 
though indirect, have two particular advantages: the measurements 
for any one source were made simultaneously, thereby avoiding 
problems with variability which can be extreme in some cases, 
and they also extended over a wide range of frequency, reaching 
as high as 43 GHz. We find our results to be consistent with the 
known counts at 30 GHz and in good agreement with the models of 
de Zotti et al. (2005) below 43 GHz; but for frequencies of 43 GHz 



and above, although there is agreement at the lower flux densities, 
our counts diverge progressively from those of de Zotti at the higher 
values, in that our predictions imply significantly fewer flat spec- 
trum sources. However, our data are sparse above a 15 GHz flux 
density of ~ 100 mJy and we caimot rule out the possibility that 
with increasing flux density there is a significant shift in the spec- 
tral index distribution towards flatter spectra, although this is not 
detected with any certainty in our measurements. 

The forthcoming wide-area survey with the Australia tele- 
scope (AT20G) should provide definitive source counts in the 
higher flux density range at a frequency of 20 GHz. Our own work 
can be seen as complementary in that it is applicable to somewhat 
lower flux densities and higher frequencies. 
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